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© A method for determining elementary circuits 
and initial values of flows in a pipe netwerk including 
one or more source nodes and a plurality of de- 
manding nodes. An upstream node and a down- 
stream node of each pipe is determined. The up- 
stream node has iess totai resistance irian the uown- 
stream node to the source node. A pipe joined 
between an upstream node and a downstream node 
is selected as an upstream pipe of the downstream 
node. A pipe which is never designated as an up- 
stream pipe is selected as a circuit origin pipe. A 
node level, which indicates the number of nodes 
included in the passage from the source node, is 
determined. An elementary circuit is determined, by 
tracing upstream nodes along a first route starting 
from a first node of the circuit origin pipe and along 
a second route starting from a second node of the 
circuit origin pipe in accordance with the node level 
of the upstream node of each pipe, until the first and 
second routes reach the same source node or the 
same demanding node. Each elementary circuit is 
used for settling initial values of flows of pipes, and 
then for computing the flows. 




FIG. 5 
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The present invention relates to a method and 
apparatus for determining elementary circuits and 
initial values of flows in a pipe network before 
applying the Hardy Cross method or the like to the 
pipe network to compute the flow of each pipe. 6 

Computation of flows in a network of pipes has 
been performed in such cases as new pipes are 
Installed or stable supply are considered in water 
or gas networks. Similar computation has also been 
carried out in electric distribution networks. io 

A method for such computation using the Har- 
dy Cross method Is described, for example, on 
pages 427 - 438 in "Fluid Mechanics". V. L, 
Streeter, at al.. McGraw-Hill Book Co.. 1985. in this 
method, an engineer must divide a pipeline net- is 
work Into a plurality of elementary circuits, and 
then designate an Initial flow for each pipe of the 
network at the Initial stage of the computation. 
Once these initial setting has been completed, the 
Hardy Cross method which has been programmed 20 
for numerical solution on a digital computer can 
accomplish the computation, and outputs the flow 
distribution of the network. Here, an elementary 
circuit refers to a loop or a pseudo loop described 
in the Streeter reference mentioned above. 25 

Determination of the elementary circuits and 
initial values, however, requires manual operation. 
This presents a problem because determining the 
elementary circuits in a complicated network of 
pipes is a very time consuming process. In addi- 30 
tion. wasteful computations are often involved in 
the process owing to engineer's errors, such as 
selecting an unnecessary elementary circuit. 

Likewise, since the initial values of flows are 
set by an Intuition of an engineer, the increasing 35 
number of iteration of computation will be required 
unless the initial setting is appropriate. This also 
presents a problem In that it takes a long time to 
complete the computation. 

It is therefore an object of the present invention 40 
to provide a method and apparatus for determining 
elementary circuits and Initial values of flows in a 
pipe network, which can automatically determine 
the elementary circuits and initial values of flows in 
a pipeline network. ^ 

In a first aspect of the present Invention, there 
is provided a method for determining elementary 
circuits and initial values of flows In a pipe network 
Including one or more source nodes and a plurality 
of demanding nodes, the source nodes and de- so 
manding nodes being interconnected with pipes, 
and each of the pipes having a known resistance, 
the method comprising the steps of: 

(a) determining a resistance of a passage from 
each of the source nodes to each one of the 5S 
demanding nodes on the basis of the resistance 
of the pipes; 




A1 2 

(b) determining, for each demanding node, a 
source node and a total resistance on the basis 
of the resistance of the passages, the total resis- 
tance being determined as the minimum resis- 
tance of the passages from the demanding node 
to the source nodes, and the source node being 
determined as a source node joined to the de- 
manding node through the passage providing 
the total resistance; 

(c) determining, for each of the pipes, an up- 
stream node and a downstream node on the 
basis of the total resistances, the upstream node 
having less total resistance than the downstream 
node; 

(d) designating a pipe joining one of the up- 
stream nodes to one of the downstream nodes 
as an upstream pipe of that downstream node; 

(e) selecting a pipe, which is never designated 
as an upstream pipe, as a circuit origin pipe; 

(f) determining, for each of the source nodes 
and the demanding nodes, a node level on the 
basis of the upstream nodes and the down- 
stream nodes by ignoring the circuit origin pipes 
of the pipeline network, the node level indicating 
the number of nodes included in the passage 
from the source node; and 

(g) determining elementary circuits, by tracing 
upstream nodes along a first route starting from 
a first node of one of the circuit origin pipes and 
along a second route starting from a second 
node of the circuit origin pipe in accordance 
with the node level of the upstream node of 
each pipe, until the first and second routes 
reach a same source node, or different source 
nodes, or a same demanding node, each of the 
elementary circuits being used as an element 
for detemnining a flow of each of the pipes. 

Here, the method may further comprise the 
steps of: 

dividing the pipeline network into the elemen- 
tary circuits by assuming that a flow of each of the 
circuit origin pipes is zero so that each of the 
elementary circuits has a tree structure; and 

determining an Initial flow of each pipe of the 
pipeline network by tracing the first route and sec- 
ond route by adding demand quantity of each node 
along the routes. 

The method may further comprise the step of 
determining the flow of each pipe in the pipeline 
network by using the initial flow of each pipe ob- 
tained in the step of determining the initial flow. 

The step of determining the flow may employ 
the Hardy Cross method. 

The pipeline network may be a gas pipeline 
network. 

The pipeline network may be a water pipeline 
network. 
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The pipeline network nnay be an electric power 
network, in which the flow is a current, and the 
pipeline is a power line. 

In a second aspect of the present invention, 
there is provided a method for determining eiennen- 
tary circuits and initial values of flows in a pipe 
network including one or nnore source nodes and a 
plurality of dennanding nodes, the source nodes 
and demanding nodes being interconnected with 
pipes, and each of the pipes having a known resis- 
tance, the method comprising the steps of: 

(a) determining a resistance of a passage from 
each of the source nodes to each one of the 
demanding nodes on the basis of the resistance 
of the pipes; 

(b) determining, for each demanding node, a 
source node and a total resistance on the basis 
of the resistance of the passages, the total resis- 
tance being determined as the minimum resis- 
tance of the ipassages from the demanding node 
to the source nodes, and the source node being 
determined as a source node joined to the de- 
manding node through the passage providing 
the total resistance; 

(c) determining, for each of the pipes, an up- 
stream node and a downstream node on the 
basis of the total resistances, the upstream node 
having less total resistance than the downstream 
node; 

(d) designating a pipe joining one of the up- 
stream nodes to one of the downstream nodes 
as an upstream pipe of that downstream node; 

(e) dividing the network into one or more tree 
structures by assuming flows of predetermined 
pipes as 0; and 

(f) determining an initial flow of each pipe of the 
pipeline network, by tracing upstream nodes 
along a first route starting from a first node of 
one of the predetermined pipes and along a 
second route starting from a second node of the 
predetermined pipe until the first and second 
routes reach the same source node, or different 
source nodes, or the same demanding node, 
and by adding demand quantity of each node 
along the routes. 

In a third aspect of the present invention, there 
is provided an apparatus for determining elemen- 
tary circuits and initial values of flows in a pipe 
network including one or more source nodes and a 
plurality of demanding nodes, the source nodes 
and demanding nodes being interconnected with 
pipes, and each of the pipes having a known resis- 
tance, the apparatus comprising: 

(a) means for determining a resistance of a 
passage from each of the source nodes to each 
one of the demanding nodes on the basis of the 
resistance of the pipes; 



(b) means for determining, for each demanding 
node, a source node and a total resistance on 
the basis of the resistance of the passages, the 
total resistance being determined as the mini- 
mum resistance of the passages from the de- 
manding node to the source nodes, and the 
source node being determined as a source node 
joined to the demanding node through the pas- 
sage providing the total resistance; 

(c) means for determining, for each of the pipes, 
an upstream node and a downstream node on 
the basis of the total resistances, the upstream 
node having less total resistance than the down- 
stream node; 

(d) means for designating a pipe joining one of 
the upstream nodes to one of the downstream 
nodes as an upstream pipe of that downstream 
node; 

(e) means for selecting a pipe, which is never 
designated as an upstream pipe, as a circuit 
origin pipe; 

(f) means for determining, for each of the source 
nodes and the demanding nodes, a node level 
on the basis of the upstream nodes and the 
downstream nodes by ignoring the circuit origin 
pipes of the pipeline network, the node level 
indicating the number of nodes included in the 
passage from the source node; and 

(g) means for determining elementary circuits, 
by tracing upstream nodes along a first route 
starting from a first node of one of the circuit 
origin pipes and along a second route starting 
from a second node of the circuit origin pipe in 
accordance with the node level of the upstream 
node of each pipe, until the first and second 
routes reach a same source node, or different 
source nodes.-or a same demanding node, each 
of the elementary circuits being used as an 
element for determining a flow of each of the 
pipes. 

The above and other objects, effects, features 
and advantages of the present invention will be- 
come more apparent from the following description 
of the embodiment thereof taken in conjunction 
with the accompanying drawings. 

Fig. 1 is a diagram illustrating a pipeline network 

to which an embodiment in accordance with the 

present invention is applied; 

Fig. 2 is a diagram illustrating a pipe Pi 3 and its 

first node FN0DE(13) and second node 

SNODE(13); 

Fig. 3 is a diagram illustrating a demand quan- 
tity q(2) of a node N2 ; 

Fig. 4 is a diagram illustrating the upstream pipe 
UPIPE(2) of the node N2 ; 

Fig. 5 is a diagram illustrating circuit origin pipes 
and node levels of the network shown in Fig. 1 ; 
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Figs. 6A - 6E are diagrams Illustrating elemen- 
tary circuits; 

Fig. 7 is a diagram illustrating computation of 
initial values of flows of the pipeline network 
shown in Fig. 1; 

Fig. 8 is a general flowchart showing stages for 
determining the elementary circuits, setting ini- 
tial values of flows, and computing the flows; 
Fig. 9 is a lilock diagram showing an apparatus 
for executing the processes of Fig. . 8; 
Fig, 10 is a flowchart illustrating the procedure 
of stage 300 for computing total resistances, 
obtaining an upstream node and downstream 
node of each pipe, and obtaining an upstream 
pipe of each node; 

Figs. 11 A and 11 B are diagrams showing a data 
table 100 containing data associated with nodes 
and a data table 110 containing data associated 
with pipes, after the first cycle of the stage 300 
has been completed; 

Figs. 12A and 12B are diagrams showing the 
data tables 100 and 110 after the second cycle 
of the stage 300 has been completed; 
Fig. 13 is a flowchart illustrating the procedure 
of stage 400 for obtaining circuit origin pipes; 
Figs, 14A. 148 and 14C are diagrams showing 
the data tables 100 and 110, and a data table 
120 containing data associated with elementary 
circuits after the stage 400 has been completed; 
Fig. 15 is a flowchart illustrating the procedure 
of stage 500 for computing a node level of each 
node; 

Figs. ISA and 16B are diagrams showing the 
data tables 100 and 110 after the first cycle of 
the stage 500 has been completed; 
Figs. 17A and 17B are diagrams showing the 
data tables 100 and 110 after the second cycle 
of the stage 500 has been completed; 
Fig. 18 is a flowchart illustrating the procedure 
of stage 600 for determining elementary circuits 
of the pipeline network; and 
Rg. 19 is a fiowchart illustrating the procedure 
of stage 700 for setting Initial values of flows of 
the pipeline network. 
The invention will now be described with refer- 
ence to the accompanying drawings. 

Before explaining the details of the present 
invention, the outline thereof will be explained with 
reference to Figs. 1 - 7. 

Fig. 1 illustrates a network of pipes. In this 
figure, nodes Ni - Nio are joined by pipes Pi - 
Pi 3. and nodes Ni and are source nodes. The 
other nodes N2, N3, N5 - Nio are demanding 
nodes. The total numbers NN and PN of the nodes 
and pipes are 10 and 13, respectively, and the total 
number SN of the source nodes is 2. The two 
source node number 1 and 4 are expressed as SN- 
(1) = 1 and SN(2) = 4. A pipe resistance RES{i) is 



given for each pipe Pj in advance. In this specifica- 
tion, the resistance of a pipe is defined as (Pi^ - 
P2^)/Q^ using the variables in Equation (1) of Chap- 
ter 2 of "Gas Engineers Handbook", C. H. Bum- 
5 ham, et al.. The Industrial Press, New York, 1966, 
Ftg. 2 illustrates a pipe Pj and two nodes at its 
ends, that is, a first node and a second node, 
which are determined arbitrary. The number of the 
first node and the number of the second node are 
10 expressed as FNODE(j) and SNODE{i). In this fig- 
ure, the pipe Pj = Pi 3. the first node number 
FNODE(13) = 9. and the second node number 
SNODE(13) = 2. The resistance of the pipe Pi 3 is 
expressed as RES(13) = 7. 
15 Rg. 3 illustrates the demand quantity q(i) of the 

node Nj. For example, the demand quantity of the 
node N2 can be expressed as q{2). 

The total resistance SD(i) from the source node 
to the node Nj is defined as follows: First, a plural- 
20 ity of passages from each one of the source nodes 
to the node Nj are specified. Second, the total 
resistance of each of the passages is calculated- 
Third, the minimum total resistance is selected as 
the total resistance SD(i). 
25 For example, the total resistance of the node 

N2 is obtained as follows: First, passages from the 
source node Ni to the node N2. and passages 
from the source node N* to the node N2 are 
specified. In this case, the pipe Pi, and the pipes 
30 Pi 2 and Pi 3 are obtained as passages from the 
source node Ni to the node N2, and the pipes P3 
and P2 is obtained as a passage from the source 
node N4 to the node N2. The total resistances of 
these passages are 10, 25, and 23, respectively. 
35 As a result, the total resistance SD (2) is 10, that is, 
the total resistance of the passage passing through 
the pipe Pi . 

Rg. 4 illustrates the upstream pipe number 
UPIPE(i) of the node Ni. The upstream pipe num- 
40 ber UPIPE(i) is specified as the number of a pipe 
. which is joined to the node Nj and which con- 
stitutes the passage associated with the total resis- 
tance SD(i) . In other words, the upstream pipe 
UPIPE(i) leads the node Nj to the source node with 
45 the minimum resistance. For example, in Rg. 4, the 
upstream pipe number of the node N2 , that is, 
UPIPE{2) = 1. 

The upstream node number UNODE{j) and the 
downstream node number DNODEO) of a pipe Pj 
50 are specified by comparing the total resistances 
SD from the source node to the two nodes at the 
ends of the pipe Pj. The node having a smaller 
total resistance is the upstream node whose num- 
ber is UNODEG), and the node having a larger total 
55 resistance is the downstream node whose number 
is DNODEQ). For example, in Rg. 1, since the 
nodes N2 and Na have total resistances of 10 and 
12. respectively, the upstream node of the pipe Pi 3 
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is the node N2, and the downstream node of the 
pipe Pi3 is the node N9. Thus. UN0DE(13) = 2 
and DN0DE(13) = 9. 

In Fig. 5. dotted lines indicate circuit origin 
pipes. A circuit origin pipe is a pipe that is never 
used as an upstream pipe. For example, the pipe 
P2 is; joined to the nodes N2 and N3, and the 
upstrean>. pipe of the node N2 is the pipe Pi and 
that of the node N3 is the pipe P3. Thus, the pipe 
P2 is not used as the upstream pipe, and hence it 
is a circuit origin pipe. Likewise, the pipes P5, P7, 
P10 and Pi 3 are circuit origin pipes. 

. . The circuit origin pipe will be described in 
more detail below. When a network of pipes in- 
cludes only one source node and no loop circuit, 
the flow of each pipe can be determined by merely 
summing demand flows of the nodes regardless of 
pressure/ On the other hand, when a network of 
pipes includes a plurality of source nodes and/or 
loop circuits^ some unknown flows arise, and hence 
it becomes necessary to determine in advance 
these unknown flows. The flows of the circuit origin 
pipes correspond to the unknown in the computa- 
tion of the pipe network. Accordingly, computation 
time will become shorter as the flows of the circuit 
origin pipes, which are assumed in advance, are 
set closer to the resultant values. For this reason, 
pipes of the least possible, flows are selected as 
circuit origin pipes, and the flows of these pipes 
are set to zero. To select the least possible flow 
pipes, the total resistance SD from the source node 
to each node is computed, and pipes joined to 
nodes whose total resistance SD is as large as 
possible are specified as having the least possible 
flows. 

After the circuit origin pipes are determined, a 
node level of each node is determined. Fig. 5 
illustrates the node level LEVEL(i) of the node Nj. 
which, is specified as the number of nodes from the 
source node to that node 'Nj. The node" level is 
specified for each ncdo in the rjct'ivorlc cf pipes 
after removing the circuit origin pipes. For exam- 
ple, the circuit origin pipes which are shown by 
dotted lines in Fig. 5 are removed from the net- 
work, and then the node levels are determined for 
all nodes as shown by numerals in small triangles. 
The node , levels of the source nodes Ni and Na 
are 1, and the node levels of the nodes N2, N5 and 
Ns next to the source node Ni are 2. The node 
levels of the nodes N3 and Ns which are joined to 
the source node N4 are also 2. and the node level 
of the node Ns which is the third downstream node 
from the source node N4 is 4. Thus, the node 
levels are determined. 

After the node levels have been determined, 
the elementary circuits are determined. First, one 
of the circuit origin pipes is selected as an origin, 
and then the nodes are traced upstream as shown 



in Figs. 6A - 6E. 

Let us consider the circuit origin pipe P2 in Fig. 
6A. The nodes at the ends of the circuit origin pipe 
P2 are nodes N2 and N3 whose levels are 2. 

5 Tracing the upstream pipes of the nodes N2 and 
Ns, the next nodes Ni and N4 whose levels are 1 
are reached. Since the nodes Ni and N4. are 
source nodes, the upstream trace is completed 
here, and one of the elementary circuits is deter- 

10 mined. That is, the circuit joining the nodes Ni . N2, 
N3 and is selected as the elementary circuit Ci 
as shown in Fig. 6 A. 

With regard to the circuit origin pipe P5 shown 
in Fig. 63, the node levels of the nodes N3 and N7 

75 at its ends are 2 and 3, respectively. The node N7 
whose node level is 3 is joined to the upstream 
node Na whose node level is 2. Then, the nodes 
N3 and Ns whose node levels are 2 are joined to 
the source node N^. Thus, the elementary circuit 

20 O2 is determined as shown in Fig. 6B. 

With regard to the circuit origin pipe P7 shown 
in Fig. 6C. the node levels of the nodes N5 and Ng 
at its ends are 2 and 4, respectively. The node Ns 
whose node level is larger than that of node Ns is 

25 joined successively to the upstream node N? and 
then to Ng whose node level is 2. After that, the 
nodes N5 and Ns whose node levels are 2 are 
joined to the source nodes Ni and N*, respec- 
tively. Thus, the elementary circuit C3 is obtained 

30 as shown in Fig. 6C. 

With regard to the circuit origin pipe Pi 0 shown 
in Fig. 6D, the node levels of the nodes Ns and 
N10 at its ends are 4 and 3, respectively. The node 
Ng whose node level is larger than that of node 

35 N10 is joined to the upstream nodes N7 whose 
node level is 3. After that, the nodes N? and N10 
whose node levels are 3 are both joined to the 
nodes N3 whose node level is 2. Since the nodes 
N7 and N10 are joined to the same node Ns, it is 

40 decided that the trace is completed here. Thus, the 

6D. 

Similar steps are taken for the circuit origin 
pipe Pi 3, and the elementary circuit C5 is obtained 

45 as shown in Fig. 6E. 

By setting the flows of the circuit origin pipes 
to zero, the network takes tree structure as shown 
in Fig. 7. In this state, the initial values of the flows 
are determined based on a demand quantity q(i) of 

50 each node. 

In Fig. 7, the nodes N2. Na, Ns. Ne, Ng and 
N10 are handled as terminal nodes after removing 
the circuit origin pipes. Therefore, the flow of each 
pipe joined to each one of the terminal nodes 

55 equals the demand quantity q of that node. For 
example, the flow of the pipe Ps joining the node 
Ns to the node N7 equals the demand quantity q(6) 
of the node Ns- The flow of the pipe Pa joining the 
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node N7 to the node Ns equals the sum q(6) + q- 
(7). Likewise, the flow of the (Dipe Pn joining the 
node Nio to the node Ns equals q(10). Therefore, 
the flow of the pipe Pg joining the node Ns to the 
source node is the sum total q(6) + q(7) + q{8) 
+ q{10). In this way, the initial values of flows of 
the pipes are set. 

After the elementary circuits Ci - Cs . and the 
initial values of flows have been determined as 
described above, the flow of each pipe of the 
network is computed by the Hardy Cross method. 
This method performs the balancing of each ele- 
mentary circuit in the network successively, and 
then repeats in an iterative manner until the sum of 
all flow corrections is less than a specified toler- 
ance. In other words, the flow of each circuit origin 
pipe is successively approximated to solutions, 
which are obtained in an iterative manner until the 
differences between the current solutions and their 
previous counterparts become less than the speci- 
fied tolerance. Details of the Hardy Cross method 
are described on pages 427 - 438 in "Ruid Me- 
chanics". V. L. Streeter. et al.. McGraw-Hill Book 
Co.. 1985. as stated before. 

The process described above will be summa- 
rized as shown in Rg. 8. 

In stage 100, known quantities are defined and 
stored in a memory. These includes the total num- 
ber NN of the nodes Ni, demand quantity q(i) of 
each node Nj. the total number PN of the pipes Pj, 
the resistance RES(j) of each pipe Pj, first and 
second nodes FNODEG) and SNODEG) of the pipe 
Pj, the total number SN of the source nodes, and 
the source node number SN(f). Here, f is 1 and 2. 
and SN(1) = 1 and SN(2) = 4 in Rg. 1, because 
the source nodes are U^ and N^. 

In stage 200, areas for unknown quantities are 
defined in a memory. The unknown quantities are 
the upstream pipe number UPIPE(i) of the node Nj. 
the total resistance SD(i) from the source node to 
the node Ni, the node level LEVEL(i), the upstream 
node number UNODEG) and the downstream node 
number DNODEQ) of the pipe Pj. the initial flow IF- 
Q) of the pipe Pj, the total number of the elemen- 
tary circuits CN. and the circuit origin pipe C(k) . 
Here, k is the sequential number of the elementary 
circuits, and takes values 1. 2, 3. 4, and 5 as 
shown in Rgs. 6A - 6E. 

In stage 300, the total resistance SD of each 
node from the source node is computed. As stated 
before, the total resistance SD{i) is defined as the 
minimum resistance in the resistances from all the 
source nodes to the node Ni. and the source node 
which gives the minimum resistance is specified as 
the source node of that node Nj. The upstream 
node number UNODEG) and downstream node 
number DNODEG) of each pipe Pj, and the up- 
stream pipe number UPIPE(i) of each node Ni are 



determined on the basis of the total resistances 
SD- 

In stage 400. the number of elementary circuits 
CN is computed. Subsequently, the pipes which 
j are not designated as the upstream pipes are se- 
lected as the circuit origin pipes. 

In stage 500, the node level LEVEL(i) from the 
source node to the node Nj is determined for each 
node Ni by using the upstream node number 
0 UNODEG) and downstream node number DNODEG) 
of each pipe Pj. 

In stage 600, the elementary circuits C(k) are 
determined by using the circuit origin pipes as 
origins and referring to the upstream nodes 
;5 UNODEG). downstream nodes DNODEG) and the 
node levels LEVEL{i). 

In stage 700. the flow of each circuit origin pipe 
is assumed to be zero so that the network is 
provided with tree stnjcture. After that, an initial 
20 value of the flow of each pipe of the network is 
determined. 

In stage 800. the flow of each pipe is com- 
puted by using the Hardy Cross method. In this 
process, the flow of each circuit origin pipe is 
25 successively approximated to solutions, which are 
obtained in an iterative manner continued until the 
differences between the current solutions and their 
previous counterparts become less than a pre- 
determined threshold. 
30 Rg. 9 is a block diagram showing an embodi- 

ment of a system in accordance with the present 
invention. 

In this figure, reference numeral 20 designates 
a CPU. The CPU 20 is connected to an in- 
35 put/output device 25. a known value storing portion 
30. an unknown value storing portion 50. and an 
application program storing portion 70. 

The known value storing portion 30 comprises 
memory areas 31 - 38 for storing the total number 
40 NN of the nodes, the demand quantity qO) of each 
node Ni. the total number PN of pipes, the resis- 
tance RESG) of each pipe Pj, a first node number 
FNODEG) of a pipe Pj. a second node number 
SNODE (j) of the pipe Pj. the total number SN of 
45 the source nodes, and the sequential number SN(0 
of the source nodes. These values have been in- 
putted in advance from the inpuVoutput portion 25 
to the areas 31 - 38. 

The unknown value storing portion 50 corn- 
so prises memory areas 51 - 58 for storing the up- 
stream pipe number UPIPE(i) of the node Nj, the 
total resistance SD(i) from the source node to the 
node Ni. the node level LEVEL(i) of the node Nj 
from the source node, the upstream node numt^r 
55 UNODEG) of the pipe Pj. the down stream node 
number DNODEG) of the pipe Pj, the initial flow IF- 
G) of the pipe Pj. the total number of the elemen- 
tary circuits CN. and the circuit origin pipe C(k). 
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The application program comprises a total re- 
sistance computing routine 71, an up- 
stream/downstream pipe determining routine 72, an 
upstream pipe determining routine 73, a circuit 
origin pipe determining routine 74, a node level 
determining routine 75. an elementary circuit deter- 
mining routine 76, an initial flow setting routine 77 
and a flow computing routine 78. 

Referring to Figs. 10-19, each of the stages 
300 - 800 will be described in more detail. 

Fig. 10 is a flowchart showing the process of 
the stage 300 executed by the CPU 20. In this 
process, data as shown in Figs. 11A and 11B are 
used. 

Fig. 11A illustrates a data table 100 associated 
with nodes Nj. The data table 100 consists of four 
columns: a column 101 for storing the node num- 
bers i; a column 102 for storing the upstream pipe 
numbers UPIPE(i); a column 103 for storing the 
total resistances SD(i) from the node Nj to the 
source node; and a column 104 for storing node 
levels LEVEL(i) of the node N;. Fig. 11B illustrates 
a data table 110 associated with the pipes Pj, The 
table 110 consists of four columns: a column 111 
for storing the pipe numbers j; a column 112 for 
storing pipe resistances RES(j); a column 113 for 
storing the upstream node numbers UNODEQ); and 
a column 114 for storing downstream node num- 
bers DNODE{i). The upstream node numbers 
UNODE(i) and downstream node numbers DNODE- 
(i) have been set in advance at the leftmost posi- 
tions in the columns 113 and 114. The data tables 
100 and 110 are set in the internal memory of the 
CPU 20. 

At step 301 of Fig. 10, each total resistance SD 
in the column 103 is set to an initial value 100 
which cannot take place actually. At step 302, the 
total resistance SD of each of the source nodes Ni 
and N4. is changed to 0.0. The source nodes Ni 
and N4. are known in advance. As a result, the 

iCiLMiv^oi iiuiiiucjio 111 uic ooiui wo \j\ mo uaia 

table 100 have been set as the total resistance SD- 
(i) as shown in Rg. 11 A. 

At step 303, switches SW1 and SW2 are set to 
0. The switch SW1 indicates that at least one of the 
initial values 100 of the total resistances still re- 
mains in the column 103. The switch SW2 in- 
dicates that the total resistance SD is changed. At 
step 304, a parameter j for selecting the pipe in 
ascending order is set at its initial value of 1 . 

At step 305. the upstream node number 
UNODEG) and downstream node number DNODEG) 
of the pipe Pj are substituted into iu and id, respec- 
tively. At step 306, the total resistances SD{iu) and 
SD (id) are compared. If the total resistance SD(iu) 
of the upstream node is greater than the total 
resistance SD(id) of the downstream node, that is, if 
SD(iu) > SD(id). the processing proceeds to the 



next step 307. whereas if SD(iu) S SD(id), it pro- 
ceeds to step 308. At step 307, the upstream node 
number UNODE(j) and the downstream node num- 
ber DNODE(j) are interchanged so that the total 

5 resistance of the upstream node becomes smaller 
than that of the downstream node. 

For example, for the pipe Pi, the upstream 
, Tibde number UNODE(I) = 1 and the downstream 
node number DNQPE(I) = 2 are obtained from the 

70 data table 110. In this case, since SD(1) = 0.0 and 
SD(2) = 100 as shown in the data table 100 of Fig. 
11 A, the total resistance of the upstream node Ni 
is smaller than that of the downstream node N2. 
that is. SD(1) ^ SD(2). As a result, the node inter- 

76 change at step 307 is not carried out. 

On the other hand, for the pipe P3. the up- 
stream node number UNODE(3) = 3 and the 
downstream node number DNODE(3) = 4 is ob- 
tained from the data table 110. In this case, since 

20 Sb(3) = 100 and SD(4) = 0.0 as shown in the 
data table 100. the total resistance of the upstream 
node N3 is greater than that of the downstream 
node N4. that is, SD(3) > SD(4). As a result, the 
node interchange at step 307 is carried out so that 

25 the node N3 becomes a downstream node and the 
node N* becomes an upstream node as shown at 
the line of the pipe P3 of the data table 110 of Fig. 
11B. 

At step 308, the total resistance SD(id) of the 
30 downstream node is compared with the initial value 
100. If SD equals 100, the processing proceeds to 
the next step 309, where the switch SW1 is turned 
to 1, whereas if SD differs from 100, the processing 
proceeds to step 310. For example, since the total 
35 resistance SD(2) of the node N2 is 1 00 as shown in 
the data table 100, the switch SW1 is turned to 1 in 
this case. At step 310, the resistance RES{j) of the 
pipe Pj is added to the total resistance SD(iu) of the 
upstream node, and the sum is substituted into 
40 SDNEW. The SDNEW is a candidate of the total 
resistariCS SD of the down stream ncde of the pipe 

Pj- 

At step 311, the SDNEW is compared with SD- 
(id). If SDNEW < SD(id). the processing proceeds 

45 to the next step 312. whereas if SDNEW ^ SD(id), it 
proceeds to step 313. At step 312, SDNEW is 
substituted into SD(id) so that SDNEW becomes a 
new total resistance of the downstream node. In 
addition, the switch SW2 is turned to 1, and j is 

50 substituted to the upstream pipe number UPIPE(id) 
in the column 102 of the data table 100. For 
example, for the pipe Pi, the resistance 10 of the 
pipe Pi is obtained from the column 112 of the 
data table 110. and is added to the total r sistance 

65 (SD = 0) of the upstream node Ni . The total 
resistance SD(2) of the node N2 is changed from 
100 to 10 in the table 100. In addition, the up- 
stream pipe number UPIPE(2) of th node N2 is 
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specified as 1 (= j) as shown in Fig. 11 A. 

At step 313. the parameter j is incremented to j 
+ 1 so that the next pipe of the data table 110 is 
selected to be processed. At step 314. j is com- 
pared with the total number PN of the pipes. If j > 
PN, that is, if all the pipes have been processed, 
the' processing proceeds to step 315. If j S PN, that 
is, if at least one pipe remains to be processed, the 
processing returns to step 305. 

At step 315. the switches SW1 and SW2 are 
tested if they are both 0. If both of them are 0, the 
process of stage 300 is completed. If at least one 
of them is not 0, the processing returns to step 
303. This is because the switch SW1 whose value 
Is 1 indicates that the total resistance SD whose 
value is 100 still remains, which means that the 
process of stage 300 has not yet completed, and 
the switch SW2 whose value is 1 indicates that the 
total resistance SD is changed, which means that it 
is not clear whether the process of stage 300 has 
been completed or not. In this case, the processing 
returns to step 303 to begin the process again. 

In the first cycle of the procedure from steps 
303 - 315, the upstream node and the downstream 
node of each of the pipes P3, Pa and Pn have 
been interchanged as shown in Fig. 11B. Further- 
more, the upstream pipe numbers UPlPE of the 
nodes N3. N7 and N10. and the total resistances 
SD of the nodes N2. N3 and N5 - N10 have been 
changed as shown in Fig. 11 A. In particular, the 
total resistance SD(3) of the node Na has been 
changed twice: one during the processing with 
regard to the pipe P2, and the other during the 
processing on the pipe P3. Similar changes are 
seen about the nodes N7 and Nio- 

Since the switches SW1 and SW2 have been 
turned to 1 at steps 309 and 312 in the first cycle, 
the processing enters the second cycle of the 
procedure from steps 303 - 315. In the second 
cycle, the upstream node and the downstream 
node of the pipe Ps have been interchanged during 
the processing of the pipe P3 as shown in Fig. 
12B. In connection with this, the upstream pipe 
UPIPE(6) of the node Ns has been changed from 
P7 to Ps , and the total resistance SD(6) has been 
changed from 35 to 23 through steps 310 - 313 as 
shown in Fig. 12A. The change in the total resis- 
tance SD was performed at step 312. and hence 
the switch SW2 was changed to 1 during this 
process. 

In this second cycle, the upstream node and 
the downstream node of the pipe Pio have also 
been interchanged as shown in Fig. 12B. In this 
case, however, no data in the data table 100 asso- 
ciated with the nodes has been changed. This is 
because SDNEW ^ SD(id) holds at step 31 1 . 

Since the switch SW2 has been turned to 1 
during the second cycle, the processing enters the 



third cycle of the procedure from steps 303 - 315. 
In the third procedure, no change in the data takes 
place, and hence both switches SWl and SW2 are 
kept 0. Thus, the data processing in stage 300 has 
5 been completed. 

In stage 300. although the upstream nodes and 
the downstream nodes are freely set at the initial 
set, they are correctly determined by the proce- 
dure describe above, leaving the results as shown 
10 in Figs. 12A and 128. 

The stage 400 will now be described below. 
The stage 400 finds pipes which do not serve as 
the upstream pipes, and sets them as the circuit 
origin pipes. 

75 Fig. 13 is a flowchart showing the process of 

stage 400 executed by the CPU 20. In this figure, 
the procedure from step 401 to 405 inverts the sign 
of the resistances RES of the upstream pipes, and 
sets the node level of each node to an initial value 
20 99, for example. The procedure from step 406 to 
step 410 selects the circuit origin pipes- 

At step 401, the node number i is set at its 
initial value 1 so that the nodes in the data table 
100 are processed in ascending order from i = 1 
25 to NN (= 10). At step 402. the total resistance SD- 
(i) of the node Nj is tested whether it is zero or not. 
if it is zero, the node Ni is a source node having no 
upstream pipe, and the processing proceeds to 
step 404. In contrast, if it is not zero, the node Nj 
30 has an upstream pipe whose number is UPIPE(i), 
and the processing proceeds to the next step 403. 
At step 403. the sign of the resistance RES of the 
upstream pipe whose number is UPIPE(i) is re- 
versed. For example, the upstream pipe of the 
35 node N2 is Pi , and hence the resistance RES of 
the pipe Pi is reversed to -10 as shown in Fig. 
148. In addition, the node level of the node Nj is 
set at 99. 

At step 404. the node number i is incremented 
40 to i + 1 so that the next node of the data table . 100 
of Rg. 14A is selected. At step 405, it is tested 
whether i > NN. If i > NN, the processing proceeds 
to the next step 406. and if i ^ NN. the processing 
returns to step 402 so as to process the next node. 
45 When the procedure from step 402 to step 405 has 
been completed, the column 104 of the data table 
100 and the column 1 12 of the data table 110 have 
been changed as shown in Figs. 14A and 148. 
Here, the value "99" is chosen because the node 
50 level can never reach 99. and so it will surely be 
changed to another value. The inversion of the sign 
of the resistance RES is for identifying upstream 
pipes, and hence it may be replaced with flags, 
each of which is attached to each one of the 
55 nodes. 

The procedure from step 406 to step 410 se- 
quentially tests the value of the resistance RES of 
each pipe Pj. If the resistance RES is negative or 
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zero, the procedure jumps to step 409. In contrast, 
if the resistance RES is positive, the procedure 
proceeds to step 408 where it increments the num- 
ber k that indicates the sequential number of pipes 
whose resistance is positive, and sets the pipe 
number j into the column 122 of the elennentary 
circuit table 120 of Fig. 14C as a circuit origin pipe. 
In this case, the pipes P2, Ps, P7, Pio and Pi 3 are 
selected as the circuit origin pipes from the data 
table 110, and are registered in the colunmn 122 of 
the data table 120. Thus, the pipes which do not 
serve as an upstream pipe are selected as the 
circuit origin pipes. 

Fig. 15 is a flowchart showing the process of 
stage 500 executed by the CPU 20. 

In this figure, the procedure from step 501 to 
step 504 sets the node level of the source nodes 
Ni and to 1. The procedure from step 505 to 
step 513 sets the node levels of the demanding 
nodes N2 , N3 , and N5 - Ni 0 - 

At step 501, a parameter f is set to 1. At step 
502, the source node SN(f) is read from the mem- 
ory, and its level is set to 1. At step 503, f is 
incremented to f + 1. At step 504. f is compared 
with SN, and returns to step 502 if f ^ SN. If f ^ 
SN, the processing proceeds to step 505. Thus, 
the node levels of all source nodes are set to 1 . 

At step 505, the pipe number j is initially set to 
1. so that the pipes in the data table 110 are 
selected in ascending order from 1 to PN. In addi- 
tion, a switch SW is set to 0. The switch SW of 
value 1 indicates the possibility that the process " 
has not yet completed, and hence the processing 
returns from step 513 to step 505. 

At step 506, the upstream node number and 
the downstream node number of the pipe Pj are 
substituted into and id- At step 507, the resis- 
tance RESO) of the pipe Pj is compared with 0.0, 
and the processing proceeds to the next step 508 
if RES{j) is negative, and proceeds to step 511 if 

origin pipe, whose resistance RES(i) has been set 
as positive or zero, it undergoes ho data process- 
ing. 

At step 508, the level of the downstream node 
whose number is id and the level of the upstream 
node whose number is ig, which are associated 
with the pipe Pj whose resistance RES{i) has been 
set as negative, are substituted to LD and LU, 
respectively. At step 508, the processing tests to 
detemnine whether the node level LD of the down- 
stream node equals 99, and the node level LU of 
the upstream node differs from 99. If the answer is 
positive, the processing proceeds to the next step 
510. whereas if the answer is negative, it proceeds 
to step 511. 

At step 510, the switch is changed to 1. and 
the resistance RESO) of the pipe Pj is returned to 



the positive value. In addition, the node level of the 
downstream node LEVEL (id) is set to LU + 1. 
Thus, the node level of the downstream node is set 
to the node level of the upstream node plus 1. For 

5 example, the node level of the node N2 is set to 2 
because the node level of the upstream node Ni is 
1. The procedure from step 506 to 512 is repeated 
for each of the pipes Pi -Pi 3 of the data table 110, 
and if j exceeds the total number PN (= 13) of the 

70 pipes, the processing proceeds from step 512 to 
step 513. 

At step 513, the switch SW is tested whether it 
is 1 or not. If SW = 1, some node Jevels may be 
undetermined. Accordingly, the processing returns 
75 to step 505, and the procedure from step 505 to 
step 513 is repeated again. 

Figs. 16A and 16B shows the state when the 
first cycle of the procedure from step 505 to 513 
has been completed. Although the node level of 
20 the node Ne . which is the downstream node of the 
pipe Ps. remains at 99. all the other node levels 
have been set. In the second cycle of the proce- 
dure, the node level of the node Ng is set to 4 as 
shown in Fig. 17A. and the resistance RES(8) of 
25 the pipe Ps is returned to positive as shown in Fig. 
17B, at step 510. In addition, the switch SW is 
changed to 1 . 

Since the switch SW = 1 when the second 
cycle has been completed, the processing enters 
30 into the third cycle. In the third cycle, the process- 
ing never passes through step 510. As a result, the 
switch SW is maintained at 0, and hence the pro- 
cedure is completed at step 514. 

After the node levels have been determined, 
35 the elementary circuits C(k) are determined in 
stage 600 as shown in Figs. 6A - 6E. 

In stage 600, upstream nodes are sequentially 
traced from a first node and a second node at the 
ends of the circuit origin pipe Pr until the routes 
40 starting from the first node and the second node 
reach the different source nodes, or the same 
upstream node including the same source node. 
The route starting from the first node at one end of 
the circuit origin pipe Pr is called an A route, and 
45 the route starting from the second node at the 
other end of the circuit origin pipe Pk is called a B 
route as shown in Rgs. 6C and 6D. Both A and B 
routes trace upstream nodes until they reach the 
same node whose node level is higher than those 
50 of the starting nodes, or the different source nodes. 

Fig. 18 is a flowchart showing the process of 
stage 600 executed by the CPU 20. 

At step 601, a memory area RNa(k, ma) for 
storing the nodes associated with the A route, and 
55 a memory area RNb(k, mb) for storing the nodes 
associated with the B route are defined. At step 
602, a parameter k is set to its initial value 1 so 
that the circuit origin pipes are processed in as- 
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cending order from 1 to LN ( = 5 in this embodi- 
ment, see, Fig. MC). 

At step 603, parameters mg and mb tor specify- 
ing the memory addresses in the memory areas 
RNa(k. ma) and RNb(k. mt) are initially set to 1. At 
step 604, circuit origin pipe Pr whose pipe number 
is COP(k) is selected, and its pipe number COP(k) 
is substituted into parameters jg and jb- For exam- 
ple, when k = 1, the pipe P2 is selected from the 
data table 120 shown in Fig, 14C, and hence ja = 
jb = 2. 

At step 605, the first node Nja and the second 
node Njb of the circuit origin pipe Pja and Pjb (ja = 
jb) are determined. In addition, the node levels La 
and Lb of the first and second nodes Nia and Nib 
are read from the data table 100 shown in Fig. 17A. 

After this, the route A and the route B are 
separately traced upstream from step 606 to step 
610. 

At step 606, the node level La of the node Nja 
of the A route and the node level Lb of the node Njb 
of the B route is compared. If the node level La is 
greater than the node level Lb. that is. if La > Lb. 
the route A is traced upstream at step 608, where- 
as if La ^ Lb. the route B is traced upstream at step 
609. Before that, at step 607, the node number ia of 
the A route and the node number ib of the B route 
are compared. If they are equal, that is, if both A 
route and B route reach the same node, the pro- 
cessing jumps to step 61 1 so that the process on 
the current circuit origin pipe Pk is completed. 
Otherwise, the processing proceeds to step 609. 

At step 608, the upstream node next to the 
node Nia and its node level are determined. To 
accomplish this, the upstream pipe number UPIPE- 
(ia) is obtained after storing the node number ia into 
the memory area RNa(k, ma), and is substituted 
into ia- Then, the upstream node number UNODE- 
Ga) of the upstream pipe Pja is determined, and is 
substituted into ia- In addition, the level of the node 
Nia IS obtained, and ma is incremented. 

At step 609, a procedure similar to that of step 
608 is executed with regard to the node of the B 
route. Thus, the next upstream node is traced for 
the A route or for the B route. 

At step 610, the node level La of the node Nia 
and the node level Lb of the node Nib are tested 
whether they are 1 . If both of them are 1 . that is. if 
both A route and B route reach one of the source 
nodes, the process for the current circuit origin 
pipe Pic is completed, and hence k is incremented 
at step 612. Othenwise, the processing returns to 
step 606 so as to trace upstream again. 

In this way, the upstream nodes are traced 
along the A route and the B route alternately until 
the two routes reach the same node (including the 
source node) or the different source nodes. With 
regard to the elementary circuits shown in Figs. 6A 



and 6C, the last nodes are the source nodes Ni 
and N4, and for the elementary circuits in Figs. 6B 
and 6E, the last node is one of the source nodes. 
With regard to the elementary circuit shown in Fig. 
5 6D, the last node is node Ns because the elemen- 
tary circuit is completed at that node. 

The procedure from step 606 to step 613 is 
repeated until all the circuit origin pipes Pk are 
processed, that is. until k exceeds the total number 
70 LN of the circuit origin pipes at step 613. 

After all the elementary circuits shown in Figs. 
6A - 6E are obtained, an initial flow of each pipe of 
the network is determined in stage 700. 

Fig. 19 is a flowchart showing the process of 
75 the stage 700 executed by the CPU 20. 

At step 701 . an area for each counted flag QF- 
(i) is defined which indicates that the demand 
quantity q(i) of the node Nj has already been count- 
ed. In addition, an area for an initial value QPO) of 
20 the flow of each pipe Pj is also defined- At step 
702. the initial values QP(j). the flags QF(i), and a 
flag FLAG are initialized to 0. and a parameter k is 
set to 1 . Here. FLAG indicates whether the proce- 
dure from step 704 to step 714 is executed with 
25 regard to the route A (FLAG = 0) or with regard to 
the route B (FLAG = 1). The procedure is ex- 
ecuted on the route A first, and then on the route 
B. The flag FLAG is used to switch the routes. 

At step 703. the circuit origin pipe COP{k) is 
30 selected, and its upstream node number ia and 
downstream node number ib are determined. After 
that, the initial flow of each pipe Pj is determined 
along the route A, and then along the route B. 

At step 704, FLAG is tested whether it is 0 or 
35 not. If it is 0, the process on the route A is started 
after the node number ia is substituted into i and 
FLAG is set to 1 at step 706. On the other hand, if 
FLAG is 1 . the process on the route B is started 
after the node number ib is substituted into i and 
40 FLAG is set to 0 at step 705. 

At step 707, a sum quantity ADQ is initialized 
to 0. At step 708. the flag QF(i) is tested if it is 1. 
which indicates that the demand quantity of the 
node Nj has been already added. If the flag QF(i) is 
45 1. the processing jumps to step 710, whereas if it 
is 0. the processing proceeds to the next step 709. 
At step 709, the demand quantity q(i) of the node 
Ni is added to the sum quantity ADQ, and the flag 
QF(i) is turned to 1 which indicates that the addi- 
60 tion of the demand quantity of the node Nj has 
been completed. 

At step 710. the upstream pipe number UPIPE- 
(i) = i of the node Ni is obtained, and the flow QP- 
(i) of the pipe Pj is added to the sum quantity ADQ. 
55 Further, the upstream node number UNODEQ) of 
the pipe Pj is set as the new node number i. Thus, 
the upstream node next to the cun^ent node Nj is 
set as the new node to be processed. 
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At Step 711, the node level of the node Nj is 
tested if it is 1 . If it is greater than 1 , the process- 
ing returns to step 708 to handle the upstream 
node and pipe. If it is 1, that is, if the node Nj is the 
source node, FLAG is tested if it is 0. If FLAG is 1. 
the processing returns to step 704 to start the 
second cycle of the procedure from step 704 to 
step 712 to perform processing on the route B. If 
FLAG = 0. k is incremented to k + 1 at step 713 so 
that the next circuit origin pipe COP(k) is pro- 
cessed; Thus, the procedure is repeated until all 
the circuit origin pipes are processed, which is 
confirmed at step 714 by testing if k exceeds CN. 

For example, for the circuit origin pipe Pio as 
shown in Fig. 7, the node Nio is selected as its 
upstream node, and its demand quantity q(10) is 
added to ADQ which is 0. Subsequently, the up- 
stream pipe Pii is selected and its flow QP{11) is 
determined by QP + ADQ = q(10)-in this case. 
Then, the upstream node Ns is selected, and its 
demand quality q(8) is added to ADQ = q(10), 
resulting in q(10) + q(a) . After this, the upstream 
pipe Ps is selected, and its flow QP(6) is set at QP 
+ ADQ = q(10) + q(8). Finally, the upstream 
node N4 of the pipe Pg is selected. 

Likewise, in the second cycle of the procedure, 
the demand quantity q(6) of the node Ng, the 
demand quantity q(7) of the node N7. and the 
quantity q(10) and q(8) are added successively. 
Thus, the initial flows as shown in Fig. 7 are 
obtained. In this case, the demand quantity of each 
node is counted only once because the counted 
flag QF(i) is set to 1 once it is added to ADQ at 
step 709. 

According to the present embodiment, the ele- 
mentary circuits necessary for computing the flows 
of a pipeline network, which includes any number 
of source nodes, demanding nodes and pipes, can 
be automatically carried out. Therefore, it can 
achieve labor saving. 

In audllion. since tlie present ernbcdirTierit sets 
initial flows of the network automatically, it can 
achieve time saving. In setting the initial values, 
which are essential for computing the flows of the 
pipeline network and have great influence on com- 
puting time, the flow of a pipe which seems to have 
the least flow in each elementary circuit is set to 0, 
so that the elementary circuit will have a tree 
structure. 

In the present invention, changes and modifica- 
tions may be made without departing from the 
scope of the invention. For example, the data in 
stage 100 can be prepared by processing and 
editing digital mapping data. If the total number of 
the nodes is not more than 50 or so, the data can 
be prepared manually. Furthermore, stage 600 may 
be incorporated into stage 800 so that the flows 
through the circuit origin pipes can be set by the 



Hardy Cross method while finding elementary cir- 
cuit by using the procedure as defined in stage 
600. 

Moreover, the present invention can be applied 
5 to the computation of a network of water pipes, of 
gas pipelines, and of an electric wires. For comput- 
ing the currents of an electric network, pipes are 
replaced with wires and flows are replaced with 
currents. 

70 The present invention has been described in 

detail with respect to an embodiment, and it will 
now be apparent from the foregoing to those skilled 
in the art that changes and modifications may be 
made without departing from the invention in its 

15 broader aspects, and it is the intention, therefore, in 
the appended claims to cover all such changes and 
modifications as fall within the true spirit of the 
invention. 

20 Clainns 

1. A method for determining elementary circuits 
and initial values of flows in a pipe network 
including one or more source nodes and a 
25 plurality of demanding nodes, said source 

nodes and demanding nodes being intercon- 
nected with pipes, and each of the pipes hav- 
ing a known resistance, said method character- 
ized by comprising the steps of: 
30 (a) determining a resistance of a passage 

from each of said source nodes to each one 
of said demanding nodes on the basis of 
the resistance of said pipes; 

(b) determining, for each demanding node, 
35 a source node and a total resistance on the 

basis of the resistance of the passages, said 
total resistance being determined as the 
minimum resistance of the passages from 
the demanding node to the source nodes, 
40 and said source node being determined as 

a source node joined to the demanding 
node through the passage providing the to- 
tal resistance; 

(c) determining, for each of said pipes, an 
45 upstream node and a downstream node on 

the basis of the total resistances, said up- 
stream node having less total resistance 
than said downstream node; 

(d) designating a pipe joining one of said 
50 upstream nodes to one of said downstream 

nodes as an upstream pipe of that down- 
stream node; 

(e) selecting a pipe, which is never des- 
ignated as an upstream pipe, as a circuit 

55 origin pipe; 

(f) determining, for each of said source 
nodes and said demanding nodes, a node 
level on the basis of said upstream nodes 
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and said downstream nodes by ignoring the 
circuit origin pipes of said pipeline network, 
said node level indicating the number of 
nodes included in the passage from the 
source node; and ^ 
(g) determining elennentary circuits, by trac- 8- 
ing upstream nodes along a first route start- 
ing from a first node of one of said circuit 
origin pipes and along a second route start- 
ing from a second node of the circuit origin io 
pipe in accordance with the node level of 
the upstream node of each pipe, until the 
first and second routes reach a same 
source node, or different source nodes, or a 
same demanding node, each of said ele- is 
mentary circuits being used as an element 
for determining a flov/ of each of said pipes. 

A method for determining elementary circuits 
and initial values of flows in a pipe network as 20 
claimed in claim 1, further characterized by 
comprising the steps of: 

dividing said pipeline network into said ele- 
mentary circuits by assuming that a flow of 
each of said circuit origin pipes is zero so that 25 
each of said elementary circuits has a tree 

structure; and 

determining an initial flow of each pipe of 
said pipeline network by tracing said first route 
and second route by adding demand quantity 30 
of each node along said routes. 

A method for determining elementary circuits 
and initial values of flows in a pipe network as 
claimed in claim 2, further characterized by 35 
comprising the step of determining the flow of 
each pipe in said pipeline network by using the 
initial flow of each pipe obtained in the step of 
determining the initial flow. 

40 

A method for determining elementary circuits 
and initial values of flows in a pipe network as 
claimed in claim 3. characterized in that said 
step of determining the flow employs the Har- 
dy Cross method. ^ 

A method for determining elementary circuits 
and initial values of flows in a pipe network as 
claimed in claim 4. characterized in that said 
pipeline network is a gas pipeline network. so 

A method for determining elementary circuits 
and initial values of flows in a pipe network as 
claimed in claim 4, characterized in that said 
pipeline network is a water pipeline network. 65 

r. A method for determining elementary circuits 
and initial values of flows in a pipe network as 



claimed in claim 4, characterized in that said 
pipeline network is an electric power network, 
said flow is a curent, and said pipeline is a 
power line. 

A method for determining elementary circuits 
and initial values of flows in a pipe network 
including one or more source nodes and a 
plurality of demanding nodes, said source 
nodes and demanding nodes being intercon- 
nected with pipes, and each of the pipes hav- 
ing a known resistance, said method character- 
ized by comprising the steps of: 

(a) determining a resistance of a passage 
from each of said source nodes to each one 
of said demanding nodes on the basis of 
the resistance of said pipes; 

(b) determining, for each demanding node, 
a source node and a total resistance on the 
basis of the resistance of the passages, said 
total resistance being determined as the 
minimum resistance of the passages from 
the demanding node to the source nodes, 
and said source node being determined as 
a source node joined to the demanding 
node through the passage providing the to- 
tal resistance; 

(c) determining, for each of said pipes, an 
upstream node and a downstream node on 
the basis of the total resistances, said up- 
stream node having less total resistance 
than said downstream node; 

(d) designating a pipe joining one of said 
upstream nodes to one of said downstream 
nodes as an upstream pipe of that down- 
stream node; 

(e) dividing said network into one or more 
tree structures by assuming flows of pre- 
determined pipes as 0; and 

(f) determining an initial flow of each pipe of 
said pipeline network, by tracing upstream 
nodes along a first route starting from a first 
node of one of said predetermined pipes 
and along a second route starting from a 
second node of the predetermined pipe un- 
til the first and second routes reach the 
same source node, or different source 
nodes, or the same demanding node, and 
by adding demand quantity of each node 
along said routes. 

9. An apparatus for •determining elementary cir- 
cuits and initial values of flows in a pipe net- 
work including one or more source nodes and 
a plurality of demanding nodes, said source 
nodes and demanding nodes being intercon- 
nected with pipes, and each of the pipes hav- 
ing a known resistance, said apparatus char- 
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acterized by comprising: 

(a) means for determining a resistance of a 
passage from each of said source nodes to 
each one of said demanding nodes on the 
basis of the resistance of said pipes: 5 

(b) means for determining, for each de- 
manding node, a source node and a total 
resistance on the basis of the resistance of 
the passages, said total resistance being 
determined as the minimum resistance of w 
the passages from the demanding node to 

the source nodes, and said source node 
being determined as a source node joined 
to the demanding node through the passage 
providing the total resistance; ,5 

(c) means for determining, for each of said 
pipes, an upstream node and a downstream 
node on the basis of the total resistances, 
said upstream node having less total resis- 
tance than said downstream node; 20 

(d) means for designating a pipe joining one 
of said upstream nodes to one of said 
downstream nodes as an upstream pipe of 
that downstream node; 

(e) means for selecting a pipe, which is 25 
never designated as an upstream pipe, as a 
circuit origin pipe; 

(f) means for determining, for each of said 
source nodes and said demanding nodes a 
node level on the basis of said upstream 30 
nodes and said downstream nodes by ig- 
noring the circuit origin pipes of said pipe- 
line network, said node level indicating the 
number of nodes included in the passage 

from the source node; and 35 
(9) means for determining elementary cir- 
cuits, by tracing upstream nodes along a 
first route starting from a first node of one of 
sajd circuit origin pipes and along a second 
route starting from a second node of the 40 
circuit origin pipe in accordance with the 
node ievei of the upstream node of each 
pipe, until the first and second routes reach 
a same source node, or different source 
nodes, or a same demanding node, each of 45 
said elementary circuits being used as an 
element for determining a flow of each of 
said pipes. 

10. An apparatus for determining elementary cir- 50 
cuits and initial values of flows in a pipe net- 
work as claimed in claim 9, further character- 
ized by comprising: 

means for dividing sad pipeline network 
into said elementary circuits by assuming that 55 
a flow of each of said circuit origin pipes is 
zero so that each of said elementary circuits 
has a tree structure; and 



means for determining an initial fiow of 
each pipe of said pipeline network by tracing 
said first route and second route by adding 
demand quantity of each node along said 
routes. 

11. An apparatus for determining elementary cir- 
cuits and initial values of flows in a pipe net- 
work as claimed in claim 10, further character- 
fzed by comprising means for determining the 
flow of each pipe in said pipeline network by 
usmg the initial flow of each pipe obtained by 
said means for determining the initial flow, 

12. An apparatus for determining elementary cir- 
cuits and initial values of flows in a pipe net- 
work as claimed in claim 11, characterized in 
that said means for determining the flow em- 
ploys the Hardy Cross method. 

13. An apparatus for determining elementary cir- 
cuits and initial values of flows in a pipe net- 
work as claimed in claim 12. characterized in 
that said pipeline network is a gas pipeline 
network. 

14. An apparatus for determining elementary cir- 
cuits and initial values of flows in a pipe net- 
work as claimed in claim 12, characterized in 
that said pipeline network is a water pipeline 
network. 

15. An apparatus for determining elementary cir- 
cuits and initial values of flows in a pipe net- 
work as claimed in claim 12, characterized in 
that said pipeline network is an electric power 

- network, said flow is a current, and said pipe- 
line is a power line. 



16. A computer program product for use with a 

SVStenn fnr riAfarry-iir^ir^^ 

...,y oionic;iiiciiy (jircuiis ano 

initial values of flows in a pipe network includ- 
ing one or more source nodes and a plurality 
of demanding nodes, said source nodes and 
demanding nodes being interconnected with 
pipes, and each of the pipes having a known 
resistance, said computer program product 
characterized by comprising: 

(a) first computer readable program code 
means for causing a computer to determine 
a resistance of a passage from each of said 
source nodes to each one of said demand- 
ing nodes on the basis of the resistance of 
said pipes; 

(b) second computer readable program 
code means for causing said computer to 
determine for each demanding node, a 
source node and a total resistance on the 
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basis of the resistance of the passages, said 
total resistance being deternnined as the 
minimum resistance of the passages from 
the demanding node to the source nodes, 
and said source node being determined as s 
a source node joined to the demanding 
node through the passage providing the to- 
tal resistance; 

(c) third computer readable program code 
means for causing said computer to deter- io 
mine for each of said pipes, an upstream 
node and a downstream node on the basis 

of the total resistances, said upstream node 
having less total resistance than said down- 
stream node; 

(d) fourth computer readable program code 
means for causing said computer to des- 
ignate a pipe joining one of said upstream 
nodes to one of said downstream nodes as 

an upstream pipe of that downstream node; 20 

(e) fifth computer readable program code 
means for causing said computer to select 
a pipe, which is never designated as an 
upstream pipe, as a circuit origin pipe; 

(f) sixth computer readable program code 25 
means for causing said computer to deter- 
mine for each of said source nodes and 

said demanding nodes, a node level on the 
basis of said upstream nodes and said 
downstream nodes by ignoring the circuit 30 
origin pipes of said pipeline network, said 
node level indicating the number of nodes 
included in the passage from the source 
node; 

(g) seventh computer readable program 35 
code means for causing said computer to 
determine elementary circuits, by tracing 
upstream nodes along a first route starting 
from a first node of one of said circuit origin 
pipes and along a second route starting 40 
from a second node of the circuit origin 
pipe in accordance with the node level of 

the upstream node of each pipe, until the 
first and second routes reach a same 
source node, or different source nodes, or a 45 
same demanding node, each of said ele- 
mentary circuits being used as an element 
for determining a flow of each of said pipes; 
and 

(h) a storing medium usable by said com- so 
puter and storing said first to seventh com- 
puter readable program code means. 



17. A computer program product for use with a 
system for determining elementary circuits and 
initial values of flows in a pipe network includ- 
ing one or more source nodes and a plurality 
of demanding nodes, said source nodes and 
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demanding nodes being interconnected with 
pipes, and each of the pipes having a known 
resistance, said computer program product 
characterized by comprising: 

(a) first computer readable program code 
means for causing a computer to detemnine 
a resistance of a passage from each of said 
source nodes to each one of said demand- 
ing nodes on the basis of the resistance of 
said pipes; 

(b) second computer readable program 
code means for causing said computer to 
determine for each demanding node, a 
source node and a total resistance on the 
basis of the resistance of the passages, said 
total resistance being determined as the 

. minimum resistance of the passages from 
the demanding node to the source nodes, 
and said source node being determined as 
a source node joined to the demanding 
node through the passage providing the to- 
tal resistance; 

(c) third computer readable program code 
means for causing said computer to deter- 
mine for each of said pipes, an upstream 
node and a downstream node on the basis 
of the total resistances, said upstream node 
having less total resistance than said down- 
stream node; 

(d) fourth computer readable program code 
means for causing said computer to des- 
ignate a pipe joining one of said upstream 
nodes to one of said downstream nodes as 
an upstream pipe of that downstream node; 

(e) fifth computer readable program code 
means for causing said computer to divide 
said network into one or more tree struc- 
tures by assuming flows of predetermined 
pipes as 0; 

(f) sixth computer readable program code 
means for causing said computer to deter- 
mine an initial flow of each pipe of said 
pipeline network, by tracing upstream nodes 
along a first route starting from a first node 
of one of said predetermined pipes and 
along a second route starting from a second 
node of the predetermined pipe until the 
first and second routes reach the same 
source node, or different source nodes, or 
the same demanding node, and by adding 
demand quantity of each node along said 
routes; and 

(g) a storing medium usable by said com- 
puter and storing said first to sixth computer 
readable program code means. 
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